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PLANETARY MIGRATION IN RESONANCE: 
THE QUESTION OF THE ECCENTRICITIES 
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Abstract. The formation of resonant planets pairs in exoplanetary sys- 
tems involves planetary migration inside the protoplanetary disc : an 
inwards migrating outer planet captures in Mean Motion Resonance an 
inner planet. During the migration of the resonant pair of planets, the 
eccentricities are expected to rise excessively, if no damping mechanism 
is applied on the inner planet. We express the required damping action 
to match the observations, and we show that the inner disk can play 
this role. This result applies for instance to the system GJ 876 : we re- 
produce the observed orbital elements through a fully hydrodynamical 
simulation of the evolution of the resonant planets. 
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CO 1 Introduction 
(N 

Z^\ Most of the known exoplanets are giant planets on short period orbits around their 
stars. In the core-accretion model, such planets should have formed much further 
I | away, beyond the so-called "snow-line" , where water is solid and contributes to the 
iy> ■ formation of a massive solid core. This suggests that inward planetary migration 
f*^ \ occured in the protoplanetary disc. Through mutual gravitationnal interactions, 
r* ■ a planet and the gaseous disc exchange energy and angular momentum. For a 
. "Z~ giant planet, this lead to the opening of a gap and to the migration of the planet 
^ ■ ( jLin fc Papaloizou, 1986"! ). 

A stronger evidence for migration is that a large fraction of multiplanct systems 
are engaged a Mean Motion Resonance (MMR) : the ratio of the orbital periods of 
the outer to the inner planet, equals that of two small integers. This cannot happen 
by chance. It betrays a convergent migration of the planets and a dissipative 
process to enter this dynamical configuration. 

For instance, in GJ 876, the two most massive planets are in a 2:1 MMR, 
with periods of « 30 and rs 60 days. To model the formation of this system, 
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Lee & Peale (2002) performed 3-body simulations of a central host star and two 
planets, with additional (dissipative) forces that reproduced the effects of disc- 
planet interaction. Only the outer planet was forced to migrate inward (assuming 
that the inner planet has no ambient gas). Capture into resonance occured when 
the outer planet crosses the location at which the mean orbital periods have a ratio 
of 2. After the resonant capture, however, the eccentricities of the two planets rise 
dramatically (except if an irrealistically high eccentricity damping is applied to 
the outer planet). 

2 Migration of a pair of giant planets : analytical calculations 

A planet of mass M p , semi-major-axis a, and eccentricity e has an energy E = 
—GM*M p /2a and an angular momentum H = M p ^/GM^a(l — e), where M* 
denote the sum of the masses of the central star and the planet. Denoting by 

r a = —a/a and r e = — e/e, one finds : 

(2.1) 
(2.2) 

If two planets are in resonance, the ratio between their energies is constant. 
We define : e — E2/E1 — MaOi / 'Mi 02 ■ The variation of the energy of the entire 
system (the pair of planets), E = E\ + E 2 = Ei(l + s), is the sum of the energy 
variations applied to both planets, such that E = E\jr ai + E2/r a2 . 

The same holds for the angular momentum. If the planets are in MMR and 
their eccentricities are constant, the ratio between their angular momenta is also 
constant : r\ = H2/H1. 

With some algebra, one finds a relation between r 0l , r Q2 , r ei , and r e2 for the 
planets to evolve in resonance with given eccentricities. It reads : 
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The damping rate that should be applied to the inner planet (l/r ei ) is the sum 
of two terms. The first one (l/r eimax ) is required to balance the action on the 
outer planet ; it is not zero in general. This explains why the eccentricities rised 
too much in Lee & Peak's simulations. 

The second term (i^i/r Ql ) is proportional to t Qi _1 , with K\ < : if the inner 
planet is pushed outwards by an inner disk (&i > 0), an additionnal damping is 
required. 

3 Inner disc 

We computed a few hydro-simulations of a giant planet, on a fixed orbit with 
various eccentricities. A gaseous disc is present inside the planetary orbit. We 
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Table 1. Parameters of the system GJ 876 as given by Laughlin et al. (2005) for GJ 876 b 



and GJ 876c, and by http://exoplanets.eu/ for GJ 876 d. 



M. sin i 
(M Jup ) 



period 

(days) 



(AU) 



1.935 
0.56 
0.023 



60.93 ±0.03 

30.38 ±0.03 

1.94 



0.20783 

0.13 
0.020807 



0.029 ±0.005 

0.218 ±0.002 

0. 



measure the otrque and power of the force of the disc on the planet. We find 
that for e > 0.1, the disc damps the eccentricity of the planet. For e < 0.3, the 
expected migration if one releases the planet is directed outwards. 

Thus, the inner disc can play a crucial role on the orbital evolution of a pair of 
planets in resonance. Note that the evolution of the inner disc depends crucially 
on the radius of its inner edge i?; n f , with respect to a, the semi-major-axis of the 



orbit of the gap-opening planet (Crida & Morbidelli, 2007) If R- ln i/a is not small 



with respect to 1, the inner disc should disappear promptly. On the other hand, 
if i?i n f <C a, it should not disappear at all. 



4 Application to GJ 876 



The two planets GJ 876b and GJ 876c are in 2:1 MMR (see table [J). A third, 
smaller planet orbits very close to the star. This planet is not massive enough 
to open a gap in the gas disc (see e.g. |Crida et al., 20 06). Thus it was in type I 
migration, and most likely migrated all the way inwards, until it falls in the empty 
cavity between the star and i?i n f. More precisely until the outermost resonance 
with the planet (the 2:1) lies in the cavity: the disc and the planet exchange 
angular momentum through the resonances. Therefore, we can say that in this 
system, i? inf was located at the 2:1 MMR with GJ 876 d, that is at 0.033 AU. 

We use a code in which the whole disc is consistently simulated down to R m t , 



thanks to the use of ID grid coupled to the classical 2D grid (Crida et al., 2007) 



We compute a simulation of the two most massive planets of GJ 876 in a disc with 
a viscosity given by an a-prescription with a = 10~ 2 and an aspect ratio of 0.07. 
First, the planets are on fixed circular orbits, and shape a gap in the disc. Then, 
they are released and they migrate freely. The results are shown on Figure [U 
A resonant capture in the 2:1 MMR happens at t = 110; then, the eccentricites 
increase, but the presence of the inner disc prevents e reaching too high values. 
In the end of the simulation, the semi-major-axes and eccentricities have almost 
exactly the observed values. 
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Fig. 1. Semi-major axis (blue, left j/-axis) and eccentricity (red, right y-axis) evolution 
of GJ 876b (light colour, oi, ei) and GJ 876c (dark, a^, e-i). The horizontal lines 
correspond to the observed values, as given by Table [T] 



5 Conclusion 

This work shows that the eccentricity of resonant giant exoplanets is governed 
by the presence of a gas disc inside the orbits of the innermost one. For given 
eccentricities, resonances and planet masses, the required action of the inner disc 
is provided analytically. In the GJ 876 system, the radius of the inner edge of 
the disc can be estimated thanks to the presence of a third planet ; a numerical 
simulation of this system reproduces very well the observed parameters, which 
confirms the validity of the model. 

In conclusion, the observed moderate eccentricities of the exoplanets in reso- 
nance suggest that in these systems, an inner disc was present during the migration. 



This work has also led to the publication of an article in A&A (Crida et al., 2008) 
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